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Abstract 

Strong - weak coupling duality in string theory allows us to compute physical quantities 
both at the weak coupling end and at the strong coupling end. Furthermore perturbative 
^ string theory can be used to compute corrections to the leading order formula at both ends. 

7^, We explore the possibility of constructing a smooth interpolating formula that agrees with 
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the perturbation expansion at both ends and leads to a fairly accurate determination of the 



»>■ quantity in consideration over the entire range of the coupling constant. We apply this to 

l/~j study the mass of the stable non-BPS state in SO (32) heterotic / type I string theory with 

(^ encouraging results. In particular our result suggests that after taking into account one loop 

■^ corrections to the mass in the heterotic and type I string theory, the interpolating function 

(y-i determines the mass within 10% accuracy over the entire range of coupling constant. 
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1 Introduction 

At present we do not have a fully non-perturbative definition of string theory except in some 
special backgrounds involving AdS spaces. As a result we can only compute perturbative 
corrections around a given background (see [I] for up to date results on superstring perturbation 
theory) which is expected to break down at finite or strong coupling. S-duality provides a way 
out at strong coupling by mapping the problem to a weak coupling problem in a dual string 
theory. However there are no known techniques for systematic computation at finite coupling. 
Certain supersymmetric quantities can be computed for all values of the coupling since once 
they are computed at weak coupling they remain valid at all couplings. This includes spectrum 
of BPS states, certain terms in the low energy effective action {e.g. the prepotential in type 
11 string theory compactified on a Calabi-Yau 3-fold) etc. Recently remarkable progress has 
been made towards determining certain other class of terms in the low energy effective action 
whose form is not protected against quantum corrections, but are sufficiently constrained by 
supersymmetry so that by knowing the perturbative answer at various ends we can completely 
fix these terms (see e.g. [2] and references therein). However this still leaves open the question of 
how to determine the wide class of other observables whose form is not in any way restricted by 
supersymmetry. Most of the interesting observables like the S-matrix and masses of non-BPS 
states fall in this category. 



In this paper we explore the possibihty that knowing the behaviour at strong and weak 
couphng and matching the results from two ends, we may be able to get fairly accurate results 
for physical quantities even at finite coupling. The idea is as follows. Let us denote by F^{g) 
and F^{g) respectively the contributions to a given physical quantity up to m-th order at the 
weak coupling end {g — )• 0) and n-th order at the strong coupling end {g — )■ oo). We then try 
to look for a smooth interpolating function Fm^n{g) whose Taylor series expansions at the weak 
and the strong coupling ends match those of the functions F^{g) and F^{g) to appropriate 
order. Under favourable circumstances the function Fm,n{g) may come reasonably close to the 
actual function F[g) for sufficiently large m, n. Since the perturbation expansion in string 
theory is an asymptotic expansion, we do not expect that we can approach arbitrarily close 
to the exact result; but the question is whether we can reach fairly close to the exact result 
following this procedure]^ 

We shall apply this procedure to study the mass of the lightest SO (32) spinor state in 
SO (32) heterotic or equivalently type I string theory. Due to charge conservation this state is 
guaranteed to be stable even though it breaks all supersymmetry. In SO (32) heterotic string 
theory this is a perturbative string state |4]-[6] whereas in type I string theory this is described 
by a stable non-BPS DO-brane [7||8] (see also p}jl2]). Thus it is meaningful to look for a 
function F{g) that will give the mass of this state as a function of the string coupling constant 

g- 

The rest of the paper is organised as follows. In [|2] we fix our conventions, describe our 
strategy for finding the interpolating function and also explicitly find the interpolating function 
at the leading order. In [|3]we compute first subleading correction to the mass of the stable non- 
BPS state in type I string theory. In ^we find the first subleading correction to the mass of 
stable non-BPS state in SO(32) heterotic string theory. In ^we find the interpolating function 
taking into account the subleading corrections at the two ends and compare the result with 
the leading order interpolating function, as well as the strong and weak coupling expansions. 
We find close matching of all these functions within about 10%, indicating that already at this 
order the interpolating function may be within 10% of the exact result for all values of the 
coupling. In appendix |X] we test the efficiency of our interpolation algorithm by applying it on 
several test functions. Appendix |B] contains some technical details of the analysis carried out 
in^ 



^In a different context but similar spirit, ref. [3] attepmted to find an approximate formula for the negative 
mode eigenvalue of the Schwarzschild black hole as a function of dimension D using the known behaviour at 
large and small £> — 3. 



2 Conventions and Strategy 



We begin by fixing the various normalization conventions we shall be using in our analysis. 
We denote by Ghiiu and Gi^^ the heterotic and type I metric, defined so that the fundamental 
strings in the respective string theories have tension 1/27? in these metrics. We shall choose 
the dilatons (pn and (pj in the two theories so that the part of the action involving the metric, 
dilaton and the SO (32) gauge fields take the form: 



Sh= f d^\e-^'t'"^/detGH 
for the heterotic string theory and 



^Rh + 2G>i,'d,^Hd,^H - ^G^j^G'^TrviF.pF, 



(2.1) 



Si= f d^\^/detGi 



-2<t>i 



Ri + 26^9^0,5,0/ !> - Ce-^^GfGf Trv(F^pF, 



(2.2) 



for the dual type I string theory. Here Try denotes trace in the vector representation of SO (32), 

C = 2-13/2^-7/2 ^2.3) 

and Rh and Rj denotes the scalar curvatures computed from the heterotic and type I metrics 
respectively. Note that the overall normalization of the terms involving the metric and the 
dilaton can be changed by shifting (pn and 0/, but once these terms have been fixed the 
normalization of the gauge field kinetic term is no longer arbitrary. For the heterotic string 
theory this normalization was determined in [6] while for type I string theory this can be found 
e.g. in |13|g 

We now introduce the Einstein metric Qf^u via the field redefinitions 



'^^Hfj.u fc- y^ivi ^1)11/ ^ 9tJ.u • 



In terms of the metric g^y the action takes the form 

Sh = / d}^x\/(ieig 
and 



\r - \g^''d,<^udAn - \-^e-^"l\^' g'^Tvy^F^.F, 



Sj= f d^^x^/d^ 



^R - -^g^^d.hdAi - Ce'^'/'g^^gP'^TrviF.pF, 



(2.4) 



(2.5) 



(2.6) 



^In 
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the action was written as a 8+1 dimensional integral in the T-dual type I' description. Here we 



have expressed the action as a 9+1 dimensional integral in type I description. 



Comparing (2.5) and (2.6) we see that the two actions agree if we make the identification 



We shall define the heterotic coupling qh and the type I coupling gj via the relations 



9H = e 



{4>h) 



9i = e 



(0i> 



Eqs.(2.7), (2.3) then give 



5^7 



gngi = 2-«C-^ = 2V 



(2.7) 



f2.8) 



(2.9) 



To test the duality between the heterotic and type I string theory, we can compare the 
fundamental heterotic string tension with the type I D-string tension. Heterotic string tension 
in Einstein frame is given by 

TH = :^igHy/\ (2.10) 



27r 



On the other hand the type I D-string tension in Einstein frame is given by 14| 15 

Td 



h2rrY/\gj)-'/\ 



(2.11) 



Note that there is an extra factor of l/-\/2 compared to that of the tension of the type IIB D- 



string to take into account the effect of orientifold projection. Using (2.9 ) we see that To = Th 



in agreement with the heterotic - type I duality 16 that identifies the fundamental heterotic 
string with the type I D-string. 

From now on all the masses will be given in the Einstein metric unless mentioned otherwise. 
In the convention described above, the mass of the lightest SO (32) spinor state in the heterotic 
string theory is given by 

2{g^Y/\l + Oigl)). (2.12) 

On the other hand the mass of the stable non-BPS DO-brane in type I string theory, transform- 
ing in the spinor representation of SO(32), can be computed as follows. To leading order in gj, 
the DO-brane is connected by marginal deformation to a D-string anti-D-string pair wrapped 
on a compact circle of radius l/-\/2 measured in the type I metric [7]. Thus its mass, measured 
in the type I metric, will be 2tc\/2 times the tension of a D-string. Converting this to the 
Einstein metric the mass of the DO-brane is given by 

2V2 7rTn{gi)-'/\l + 0{gj)) = 2' n'/' {gi)-'/\l + 0{gj)) = 2~'/'n~'/' {gH^'il + 0{gjl)) . 

(2.13) 



Comparing ( |2.13 ) and (2.12) we see that the leading order heterotic and type I results meet 
at gn = (27r)^/^. In view of this we introduce a rescaled coupling parameter g via 



g^ = 2^^/^ g, = 2'/V/'g-' 



(2.14) 



so that the two formula meet at g = 1. Furthermore since at the meeting point the mass is 
given by 2^^/^7r'^/^, we defined a renormalized mass function F{g) via the relation 



M{g) = 2''^V/'Fig) . 



(2.15) 



In terms of g the leading order weak and the strong coupling formulae (2.12) and (2.13) for the 
renormalized mass function F{g) can be expressed as 



Fo^{9)=9 



1/4 



Fo'ig) = 9 



3/4 



(2.16) 



Our goal will be to explore to what extent we can determine the full function F{g) by finding 
an interpolating function that matches onto the above functions (and perturbative corrections 
to them) at the two ends. 

Let us denote by F^{g) and F,^{g) the formulae for F{g) in the weak and strong coupling 
limits to m-th and n-th order in expansion in powers of g and g~^ respectively. In that case 
we shall choose our interpolating function aqj 

l/{2(m+n+l)} 



Fm,ni9) = 9 



1/4 



l + a,g + ---amg"' + br,g'^+' + h 



^1^+' + 



+ 6i^'"+" + ^'"+'^+1 



(2.17) 

By construction this formula reduces to the correct forms given in (2.12) and (2.13) in the 
weak and strong coupling limits. The coefficients ai, • • • a^ are determined by demanding that 
the expansion of -Fm,n in powers of g matches the weak coupling perturbation expansion to 
m-th order, while the coefficients bi,- ■ -bn are determined by demanding that the expansion 
of Fm,n in powers of 1/g matches the strong coupling perturbation expansion to n-th order n 
Since in the heterotic string theory the expansion is actually in powers of g"^, this will imply 
vanishing of the a^s for odd m. 

Clearly many other interpolations are possible. In particular we could use (fractional) power of a rational 
function for this purpose. In each case we need to determine the efficiency of the interpolating algorithm by 
studying its convergence properties. As we shall see, for the problem of studying the mass of lightest SO (32) 
spinor states in heterotic string theory, (2.171 seems to give reasonable results. 



"^Note that the a^'s and 5fc's which appear m the analysis of each F,n,n are different and have to be determined 
afresh every time. 




Figure 1: Graph of tan"^ F{g) vs. ta.n^^ g for F = Fq (thin sohd curve), F^ (thin dashed 
curve) and the interpolating functions Fqq (thick dashed curve) and Fi o (thick sohd curve). 
This should be compared with Fig. |2] which describes similar curves after taking into account 
the first subleading corrections at both ends. 



Note that this procedure is not foolproof since the term inside the square bracket in (2.17) 



could become negative and hence the right hand side of (2.17) will cease to give a real function. 



This will signal breakdown of this procedure. However for sufficiently smooth functions we 
could expect such negative a^ and/or b^ coefficients to be small even if they are present. In 



this case the term inside the square bracket in (2.17) will remain positive for all positive g and 



the procedure should continue to work. Nevertheless this approach will clearly be insensitive 
to terms in F(g) whose Taylor series expansion vanishes at both ends, e.g. e~^^~^^^ times 
any polynomial in g,g~^ for positive constants A, B. More generally since the perturbation 
expansion in string theory is expected to represent asymptotic series at both ends, we do not 
expect to get arbitrarily close to the exact result by going to arbitrarily high order. Typically for 
any given value of g the best result will be reached at some particular order in the perturbation 
theory. The hope is that this approach may take us sufficiently close to the exact formula over 
the entire range of g. In appendix |A] we have tested this procedure on several test functions. 



For example for the leading order functions (2.16) we have the interpolating function 

F,M=9"\l + 9Y". (2.18) 

Also knowing that the weak coupling expansion begins at order g"^ so that ai = 0, we get 



F,M=9''\^ + 9 



2N1/4 



(2.19) 



In Fig. LLlwe have plotted the functions F^{g), F^ {g) = F{^{g) as well as FiQ{g) and -fo,o(5') 



3 Strong coupling expansion 

We shall now compute Ff{g) by determining the first order correction AM to the mass formula 
from the strong coupling end, i.e. in type I string theory. This is given by open string one loop 
correction to the energy of the non-BPS DO-brane of type I string theory, and can be expressed 

as fnpi 



- AM = - g]'^ (Stt^)"^/^ / s~^/^ ds [Zns^dodq - Zr-dodo + Zns;D0D9 - Zr.dod9] (3.1) 
^ Jo 

where Z^s;Dodq, Zr-dq^q, Z^s-^dod^, Zr.dqdq denote respectively the contributions from the 

NS and R sector open strings with both ends on the DO-brane and NS and R sector open 

strings with one end on the DO-brane and the other end on the D9-brane. Explicit calculation 

gives 17 19 



^fiiqY f2{iqYfM) f2{iqfh{iq) 

lf2{qf 



Zr-dodo 



^NS;D0D9 



^R;D0D9 



2/i(gT' 

^Q^fiiq)f2{q? 



16v/2 



fM'fsiq) ' 
fi{q)f3{qr 
Uqff2{q) ' 

(3.2) 



where 



q = e-^' , (3.3) 

oo 

fiiq) = q'/''l[{l-q'n=vi'2r), g = e^^^ 

n=l 

oo 

n=l 

oo 

Mq) ^ q~'^''U(^ + q'''-') = v{2rr/{v{^rUr)), 



n=l 



^The variable s is related to the variable t appearing in eq.(9) of il4j a.s t — 47rs. There is a factor of 1/2 
in our expression compared to that in 14^ since we are considering the self energy of the DO-brane, and hence 
we do not get a factor of 2 in the spectrum by exchanging the two ends of the string. The explicit form of the 
open string partition function is of course different since we have a non-BPS DO-brane in type I string theory 



while 14 was considering BPS D-p-branes in type II string theory. 



U{q) 



= ^-1/24 



n(i- 



,2n-l^ 



r/(r)/r/(2r) 



n=l 



(3.4) 



Individual terms in (3.1) are both infrared (IR) and ultraviolet (UV) divergent; so we need 
to be careful about using the right prescription for the IR and UV regularization. Since 2tts 
denotes the proper length of the open string progagator flowing in the loop in the Schwinger 
representation, the IR regularization is done by putting a uniform upper cut-off A on s in 
all the integrals. Regulating the UV divergence is more subtle. As we shall describe shortly, 
instead of putting a uniform lower cut-off on the integrals, the correct prescription is to use 
a lower cut-off e on all the integrals involving Zr-dodo, Zms,dod9, Zr-dqdq and the first term 
in Zi\is;DODo, but use a cut-off e/4 on the integrals involving the second and third terms of 
Zns;Dodo- Thus the integral appears as 



-AM = -K,{gj) 



1/4 



(3.5) 



where 



K. 



- (87r^)"^/2 lim lim 

2 A^oo e-^0 



\-3/^..ri/3(.1 



1 /2(97 

2/i(g)« 2h{qY 



+ IQV2 



;^\9 



fi{q)f2{q) 
UqYUq) 



A 



'e/4 *- 



fi{qff2{q) J 
f2{iqyUiq)' 



2^/2(1 + 



hiiqfhiiq) \ 

f2{iqfh{iq)^ 



(3.6) 



Using known modular transformation laws of the /j's we can also express (3.6) in the 'closed 
string channel': 



i^. 



where 



lim lim— (87r^)"^/2 

A^oo e^O 47r 



Coo 
M 

C'o9 

q 



dt{Coo + Co, + C*,) + I dt{M 

7r/A Jir/iA 



M*] 



(3.7) 



7r\9/2 

1) 



fsiqf h{qf 



29/2 

29/2 



Ui{qf h{qf\ ' 
h{nffi{iq) fSqffSq) 

f2{^qfU^q) f2{^q?h{^q) 
UqffM UqffM ' 
Vf2{qffM f2{qYh{q)\ ' 



= e 



(3.^ 



In going from (3.6) to (3.7) we have made a change of variables t = n/s in the first integral 
and t = IT /4s in the second integral, explaining the limits of integration in (3.7). Physically 



Coo denotes the cylinder amplitude with both boundaries lying on the DO-brane, given by 
the inner product between the boundary states of the DO-brane, Ai denotes the Mobius strip 
amplitude with the boundary lying on the DO-brane, given by the inner product between 
the boundary states of the DO-brane and the crosscap, and Cog is the cylinder amplitude 
with one boundary on the DO-brane and the other boundary on the D9-brane, given by the 



inner product between the boundary states of the DO-brane and the D9-brane. In (3.7) the 



parameter t has the interpretation of the proper length of the closed string propagator in the 
Schwinger representation, and hence the upper cut-off on t should be a uniform number for 



all the integrals. This is indeed the case in (3.7) since all the upper cut-offs are vr/e, but this 
required choosing the lower cut-off on the s integral in precisely the way we have chosen in 



(3.6) 



Expressing (3.5) in terms of g defined in (2.14), carrying out the rescaling given in (2.15) 



and combining this with the leading order strong coupling result given in (2.16) we see that 
the strong coupling result for F{g) corrected to first order in 1/g is given by 



Ff{9)=9'^Hl + Ks 



K. = 2-^/^K. . 



(3.9) 



Numerical evaluation of the integrals appearing in ( |3.6 ) gives 

K, ~ .351 . 



(3.10) 

A graph of tan~^ ^lis) vs. tan^^ g can be found in Fig. M but we shall postpone the analysis 
of this function till ^ by which time we shall also compute the first subleading corrections 
from the weak coupling end. 



4 Weak coupling expansion 

Next we shall compute the first order correction to M in the weakly coupled heterotic string 
theory leading to an expression for F^{g) = F^{g). The left-handed part of the state is 
created by the spin field of SO (32) acting on the vacuum. The right handed part can be 
chosen in different ways. In the light cone gauge Neveu-Schwarz-Ramond (NSR) formalism 
the NS sector states can be taken to be of the form [5] 



C3/2I0), ^: 



-l/2V'ii/2V'-l/2|0), 



^^/2"ii|0) 



(4.1) 



10 



where '?/'L„ and aL„ are transverse fermionic and bosonic oscillators respectively with 1 < i < 8. 
This gives a total of 8+56+64=128 states. Together with the 128 fermionic states arising from 
the R-sector, these form a single long supermultiplet transforming in the spinor representation 
of the gauge group SO (32). Thus all the states suffer the same mass renormalization and we 
can focus on any one of them for computing the correction to M. We shall use covariant NSR 



formulation 20 for our computation and choose the unintegrated vertex operator in the —1 
picture to be 

icce-''S^{z)i;\z)i;\z)ij'^{z)e"''>^" (4.2) 



where (p is the bosonized ghost of the /3-7 system 20 , b, b, c, c are the diffeomorphism ghosts, 
Sa are the spin fields of SO (32) and k^ is the energy of the particle in the rest frame given by 
M. The unintegrated vertex operator in the zero picture is then given by 

ccVoik, z) = icce"'''^' S^iz) U\z)ij^iz)dX%z) + i;\z)ij'^iz)dX\z) 

+^P''{z)^\z)^X\z) + zkoij\z)ij\z)ij^{z)ij''{z)] , (4.3) 



plus terms of order e'^ whose correlation functions vanish by if-chaige conservation. For com- 
puting the torus amplitude we need to convert one of the vertex operators to integrated vertex 
operator Vq by removing the cc factor. If (5M^ denotes the one loop correction to M^, then up 
to an overall multiplicative factor SM"^ is given by 



6M' ~ {gHY / d'r / d'z {Voi-k, 0)Voik, z)) matter Z ghost , (4.4) 

where r is the modular parameter of the torus, integrated over the fundamental domain, 
( ) matter dcuotcs Correlation function in the matter sector on the torus multiplied by the matter 
partition function and Zg^ost denotes the ghost partition function after removal of ghost zero 
modes. 

To evaluate this contribution, let us first focus on the contribution from the holomor- 
phic fermions ip^^{z). The possibly non- vanishing correlation functions are of two types: 

(V'^(0)?/;2(0)^i(z)^2(^)^an^j^gpgrniutationsand(^i(0)^2(^0)V'^(0)V'°(0)V'H^)V'^(^)V'^(^)V'°(^))- 
To evaluate these we first perform a double Wick rotation on ip^ to make it '?/'^(0) correspond- 
ing to some Euclidean direction 4, at the cost of picking up a factor of i. Next we introduce 
complex fermions x^ via x^ = (^'^ + i ip'^^'^) / \/^ so that we can express tp^ as (x^ + x'^)/v^ 
for 1 < A; < 4. Then we have 

(^i(0)^^(0)^i(z)^2(z)) = -{^\Q)i,\z)i,\Q)i,\z)) 

11 



+ {x\0)x\z)fiO)x'iz)) + {x\0)x\z)fiO)x'iz))) . (4.5) 

Now for a single complex ferinion like x^{z) the correlation function on the torus is given by, 
up to a phase, 

{X\0)x\z)) = ivir))-' -ffiM^), {X\0)x\z)) = iv{r))-' -fjiM-^) , (4-6) 

Vli[Z) vu[Z) 

where z/ denotes the spin structure on the torus taking values 00, 01, 10 and 11 and t?,^ are the 
Jacobi theta functions. After combining the contribution from all the holomorphic fermions 
and the superconformal ghosts, and taking into account the r-dependent normalization factors 
the result is, up to an overall constant factor, 

{i;\0)i;'{0)i;\z)i;'iz)) 

(4.7) 
where 6,^ = 1 for z/ = 11 and 00 and 6,y = —1 for z/ = 10 and 01. The extra factor of 1/2 in 



(4.7) comes from the GSO projection in the right-moving (holomorphic) sector leading to the 



sum over spin structures u. Using the Riemann identity 

^5^79^(2;i)^^(z2)^4-23)^^(-24) = 2^ii((zi + Z2 + Z3 + Z4)/2)i?n{{zi + Z2 - Z3 - Z4)/2) 

V 

^11 ((2:1 - 2:2 - ^3 + 2;4)/2)^ii((2i - ^2 + ^3 - 2:4)72) , 

(4.^ 



and that 'd\\{^^ = we see that the right hand side of (4.7) vanishes. 



Thus we are left with the two point function of the operators appearing in the last term 



inside the curly bracket in (4.3) 



(A;o)'(e-^^°^"(°)e^'=°^"(^)^"(0)^,(^)V^i(0)V^2(0)^3(0)^^(0)V^i(^)V^2(^)^3(z)^^(z))™„,ie.^,hos* 

(4.9) 
where we have taken into account an extra — sign from the Wick rotation taking tl)^ to ?^^. Now 
following the same method as described above, the net contribution from the ten holomorphic 



12 



fermions and the superconformal ghosts to (4.9) is given by 
1 



24 



(.(r))-'""^'^^^ 



^ii(^) 



J2 Su (Mz)' + AM^fM-^) + QMzfA 



-z 



+AMz)M-zf + M-zT)- 



(4.10) 



We could manipulate this further using the Riemann identity, but for reasons that will become 
clear in appendix [B] we shall postpone this till the end. 

The rest of the contribution can also be evaluated using standard method. Using the known 



correlator between the spin fields 21 we get the contribution from the 32 anti-holomorphic 
fermions to be 



1 



^n(O) 



^ V *'ii(z) 



Y.^A^m 



16 



(4.11) 



where the factor of 1/2 now arises from the GSO projection on the left-handed (anti-holomorphic) 
fermions. Finally the contribution from the ten scalars corresponding to the space-time coor- 
dinates and the diffeomorphism ghosts together give, up to a normalization. 



^2 



W) 



16 



exp[-47r2;2/r2] 



^ii(^) 



^'ii(O) 



(4.12) 



where we have used the on-shell condition (/cq)^ = 4. Substituting (4.10), (4.11), (4.12) into 



(4.9), and using the known relation ^n(O) = — 27rr7(r)'^, we can now write down the general 
formula for 5M^: 

5M2 = 1 AT {guf {kof f rfV fd'z[{ (r^(r))-4^n(^)"^ J] M^W} 

v' 

\Y.^- (^- W + ^^v{z)H,{-z) + m,{z)H,{-zf + 4^,(z)t9,(-z)3 + d,{-zf^ 

V 

((r^(r))«(^n(^))-'}{(^(r))-^^(^)-^^(^n(^))^(M^)^exp[-47r^2/^2](r2)-^ 



(4.13) 

where M is an overall normalization constant. In this expression the factor inside the first curly 
bracket gives the contribution from the 32 left-moving fermions, the factor inside the second 
curly bracket gives the contribution from the 10 right-moving fermions and the commuting 
superconformal ghosts, and the factor inside the third curly bracket gives the contribution 



13 



from the ten scalars and the diffeomorphism ghosts. The normahzation factor A^ has been 
computed in appendix [B] by comparing the result with the expected result in the effective field 
theory. The result is 



M 



1 



210^^8 



(4.14) 



Using the fact that all the dy except du are even under z — t- —z, and dii{z) is odd under 



z — )■ —z we can express the term in the second line of (4.13) as 



16 5^ M. 



16^ii(^) 



■16^n(#, 



(4.15) 



where in the last step we have again made use of the Riemann identity (4.8). Finally using 
the fact that k"^ = M^ to leading order, and using (2.14), we get 



6M = MK^g\ 



(4.16) 



where 



K., 



~ j d\ j <fzY,[dAz/2Y^]Hr)V\v{r)r 






exp[-47rz2/T-2] (7-2) 



(4.17) 

It is easy to verify that the integrand is invariant under 2;— j-z + l, z -^ z + t and r — t- r + 1 
and (z, r) — )■ {z/t, — 1/r). The domain of integration over z ^ Zi + iz^ and r = ri + iT2 is the 
fundamental region which can be taken to be < 2:1 < 1, < 2:2 < '^"2) and the regin in r plane 
bounded by the curves t\ = ±1/2 and |r| = 1. Numerical evaluation of the integral gives 



K ~ 



23 



(4.18) 



Using (2.14), (2.15 ), the weak coupling expansion of the mass formula F{g) up to order g^ now 

F^ig) = FY{g) = g'l' (i + K^g^) ■ (4.19) 



takes the form 



5 Analysis of the results and interpolating function 

Let us first summarize the results of ^ and ^ We have found that the weak and strong 
coupling expansions of the mass function Fi^g^ are given by 

F^{g)=g^l\\^K^g'^0{g'y), F^{g) = g''\l+K,g-^+0{g-^), K^ c^ .23, K, c^ .351 . 

(5.1) 
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Figure 2: Graph of tan ^ F{g) vs. tan ^ g for F = Ff (thin sohd curve), F 
(thin dashed curve) and the interpolating function F = F3 i (the thick sohd curve). 



Ff( 



FD 




Figure 3: Graph of F2 (g) / F3^i{g) (dashed curve) and F^ (g) / F3^i{g) (continuous curve) vs. 
tan~^ g. 
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1.20 




1.15 
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1.10 


^ ^ ~- 




..'■ / \ 


1.05 


- ,■■•■'/ 


LOO 
0.95 


^.^-^"^^^ ' ""^^^'-^-C:^ 




0.90 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



Figure 4: Graph of F^^J^g) j Fj^^xi^g) vs. tan~^ (7 for various [171,11). The labels are as follows: 
thin dots for Fqq, thick dots for Fi 0, small thin dashes for ^2^0? small thick dashes for F3 0, large 
thin dashes for Fq^i, large thick dashes for Fi^i, continuous thin line for F2^i and continuous 
thick line for Fsi. 



Note that both Kg and Ky^ are smaller than unity. This implies that the corrections from both 
ends remain smaller than unity at g = 1 where the leading order weak and the strong coupling 
curves meet. Given this it is not unreasonable to expect that string perturbation theory may 
be able to give a fairly good result for the function F{g) over the entire range of g. 

We shall now follow the procedure outlined in ^ to find the interpolating functions whose 
Taylor series expansion around g = and/or g = 00 agree with those of F^{g) and/or F^{g). 
We find 

FoM=g'^\l + AKsg + g'Y/\ (5.2) 

FxM=9'^H^ + GK,g' + 9'Y^\ (5-3) 

F,M=9'^Hi + QK^9' + 9'Y^\ (5.4) 

F,M = 9"' (1 + 8 K^g" + Sir,/ + g'^fl'' , (5.5) 



i^3,o(/ =//'(! + 8 ir„(72 + ^ 



4M/ 



and 



F^M=9'" (l + 10ir^/ + 10ir,(7^ + /) 



5\1/10 



(5.6) 



(5.7) 



Fig. shows the weak and strong coupling results F^ and Ff including the first subleading 
corrections and the interpolating function F3_i((7). FigJ3^ shows the ratios F^ {g) j F'i^x{9) and 
F^ {g) / F^^i{g) as a function of tan~^ (7. We see that these ratios remain close to unity during 
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(m,n) 


(0,0) 


(0,1) 


(1,0) 


(1,1) 


(2,0) 


(2,1) 


(3,0) 


(3,1) 


F^,„(l)/F„,„(l) 


-0.125 


-0.103 





-0.055 


-0.017 


-0.066 


0.010 


-0.006 



Table 1: Table showing the value of F^„(l)/-Fm,n(l) for different {m^n). 



most of the range of g except at very large and small values of g where F^ {g) and Ff{g) 
respectively are clearly bad approximations to the actual function F{g). More specifically, if 
we denote by gc the point where F^{g) and Ff{g) meet, then F^{g) agrees with F^ilg) below 
gc within 8% and Ff{g) agrees with F^^i{g) above gc within 8%. 

To estimate how close we are to the best possible approximation to the full function F{g), 
we can analyze how the interpolating functions at different orders differ from each other. To 
test this we have plotted in Fig. Ill the ratio of Fm,n{g) / F^,i{g) for 0<??i<3, 0<r2<lasa 
function of tan^^ g. As we can see, most of the ratios remain within 10% of unity throughout the 
whole range oi g, indicating that F^^i{g) may be within 10% of the best possible approximation. 
In particular we see that F^i and F2^i lie within 5% of each other over the entire range of g. 

Another crude test for determining how good the interpolation formula is its smoothness. 
This in turn can be determined by computing F'^ ^{g) / Fm,n{g) around the matching region 
(7 ~ 1. In Table fll we have shown the value of F^^„(l)/Fm,n(l) for various interpolating 
functions. As we can see, this ratio is largest for the functions Fq^ and Fo,i - precisely the two 
functions whose deviation from F^^i is maximum in Fig. Ill This gives us indirect indication 
that -Fm.n for {fn^n) other than (0,0) and (0,1) are smoother, and hence are likely to be better 
approximations to the actual result than Fq^q or Fq^i. 

Acknowledgement: I would like to thank Rajesh Gopakumar, Dileep Jatkar and Barak 
Kol for discussions and their comments on the original manuscript. This work was supported 
in part by the J. C Bose fellowship of the Department of Science and Technology, India and 
the project 11-R&D-HRI-5.02-0304. 



A Testing the algorithm with test functions 

In this appendix we shall test the efficiency of the algorithm outlined in [|2] by applying it on 
some test functions. 
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• The first function we consider is 

F{g) = {l + g + g'fl\ (A.l) 

Tlie successive approximations Fm,n{g) to tlie function F{g) are taken to be of the form 

{l + aig + --- a^g"" + M'"+^ + • • • fei^^^" + ^™+«+i)i/(™+"+i) ^ (A.2) 

where the coefficients ai,---am are fixed by requiring that the Taylor series expansion of 
Fm,n{g) around g = matches those of F{g) up to order g"^, and the coefficients 61, ■ ■ ■ 6„ are 
fixed by requiring that the Taylor series expansion of Fm^n{g) around g = 00 matches those 
of F{g) up to order g^^''^^. For m + n + 1 even, Fm,n and F are exactly equal since we get 
l + aig + --- a„^™ + M™^^ + • • • fei^^^" + ^'"+"+^ = (1 + ^ + ^2)(™+n+i)/2, p^^ odd m + n + 1 
we get 

FoM = 1 + 9, F,M = (1 + 3g/2 + 3^/2 + g^/^ 

F2A9) = (1 + %/2 + 35(?V8 + 35//8 + 5(?V2 + /)'/', 

F3,3(^) = (1 + 7g/2 + 63g^/8 + 175^716 + 175//16 + 63^78 + 7//2 + ^^)^/^ 

F^^^{g) = (1 + 9^/2 + 99^78 + 357^716 + 3843^/128 + 3843^/128 + 357//16 



+99^/8 + 9//2 + ^9)^/^ 
etc. Using these we find that over the entire range of g 



(A.3) 



-TQ.O 

F 

-^3,3 



< .155, 

< .00043, 



^1, 



4,4 



< .013, 



< .0001, 



2,2 



< .0021, 



(A.4) 



etc. Thus the error decreseas for large m, n. 
• The second test function we consider is 



F{g) = 2g^/\l - e-^'^f'\l + e-^Y^ + g'^'^e'^'^ 



(A.5) 



Note that there are non-perturbative corrections at both ends and hence we do not expect to 
approach the exact result even by going to arbitrary high order. The approximation -Fm,n(fi') 
is again taken to be of the form 



g^'*{l + aig + ---amg"^ + hng 



m+l 



t m+n _|_ m+n+l\l/2{m+n+l) 



(A.6) 



^0,1 (^) 

FiM 
F2M 

FsAa) 
F^g) 



The successive approximations are given by 

i^o,o(^) = g"\i + gY" 

= g'^\i + ^g + g'Y^' 

= g'/\l + 3g + 6g' + g'Y/' 

= g^/\l + 5g + 45^74 + 45^''^ + 10/ + g''^^^^ 

= ^^/^(l + 7g + 91^74 + 539//I2 + 2905//8 + 91^^ 

= ^^/^(l + 9^ + 153^74 + 405//4 + 1479//8 



Ug^' + g 



7^1/14 



+97323//32 + 6519//8 + 153/ + 18/ + g 



9\1/18 



(A.7) 



etc. From this we find that over the entire range of g 



Fo,o 
F 


-1 


< .119, 


^1,1 
F 


-1 


< .064, 


F2,l 

F 


-1 


< .077, 


-^2,2 
F 


-1 


F 


-1 


< .043, 


F 


-1 


< .057 













< .037, 



(A.8) 

etc. This shows that while the F^^nS can take us quite close to the actual result, we do not 

approach arbitrarily close to the exact result by going to higher orders. The best result is 

obtained for -^2,2(5') which comes within 4% of F{g) over the entire range of g. 

• So far the functions we have analyzed have all the aj's and 6j's positive. We shall now give 

the example of a function that has some of these coefficients negative but small enough so that 

the approximation scheme used here still takes us sufficiently close to the exact function. We 

consider 

F(/ = (l-/5 + /)V2. (A.9) 

The -Fm.n's are taken to be 

(1 + ai^ + ■ ■ ■ amg"" + bng"'^^ + ■■■ 61^™+" + ^™+"+i)i/(-+n+i) . (A.io) 

In this case again for m + n + 1 even, Fm,n and F are exactly equal. For odd ?7i + n + 1 we get 

Fo,o(/ = l + g, F,,i{g) = (1 - 0.3g - 0.3/ + /)l/^ 

F2,2(^) = (1 - 0.5^ + 2.575/ + 2.575/ - 0.5/ + /)^/^ 

F3,3(^) = (1 - 0.7^ + 3.675/ - 1.7675/ - 1.7675/ + 3.675/ - 0.7/ + g^^^, 
F^^^{g) = (1 - 0.9^ + 4.815/ - 3.2025/ + 8.66644/ + 8.66644^^ - 3.2025/ 

+4.815/- 0.9/ + /)^/^ (A.ll) 
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etc. Using these we find that over the entire range of g 



Fpfl 
F 

-^3,3 



-1 
-1 



<.5, 



< .08, 



^1, 



F 

Fa A 



- 1 

- 1 



<.17, 



2,2 



< .08, 



< .05, 



(A.12) 



etc. However the error does not reduce uniformly; instead it fluctuates around the 5% mark 
as we go to higher order. Since the sign of the error fluctuates, we can considerably reduce the 
size of the error by averaging the result for different values of m, n. 
• Finally we consider the function 

/oo 
dxe-"'/'-^'"'. (A.13) 

-oo 

The expansion around (? = is known to be asymptotic. Since the function goes as A/27r as 
(? — )• and as (45f)^^/^r(l/4) as (7 — )■ oo we take 

>, -l/{2(m+n+l)} 
+ ^8^)(m+n+l)p(^/4)-2(m+n+l)^m+n+ll ^ (^^^4) 

and adjust the coefficients a^ and hk so as to reproduce the series expansion around g = ^ and 
g = oo. Note that in this case the weak coupling expansion is in powers of g^ whereas the 
strong coupling expansion is in powers of 1/g, exactly as in our case. We find that Fm,n{g) 
approaches quite close to F{g) for large enough [171,71). For example 



F4,4(^) = V2^\ 1+54^+594^+780. 788 /+1294.34/+1475.35/+1038.59/+341.428c/^ 



-1/18 



(A.15) 



agrees with F{g) within 0.4% over the entire range of g. 



B Determination of the normalization constant Af 



In this appendix we shall determine the normalization constant A/" appearing in (4.13) by 
comparing this with the corresponding expression for 6M'^ in the field theory limit. Let Xs 
denote a massive scalar in the spinor representation of the SO (32) gauge group. Then the 
action of the low energy effective field theory describing the SO (32) gauge fields and their 
coupling to the scalars Xs is given by 

1 



d'\ 



n9H) 



,Trv{F,,F^n - D.xlD^Xs - M'xlx 



(B.l) 



20 



\ 



Figure 5: A Feynman diagram of low energy effective field theory contributing to mass renor- 
malization. The thick line represents the propagator of the massive scalar in the spinor repre- 
sentation of SO (32) and the thin line represents a gauge boson propagator. 



Note that the gauge field action has been normalized in accordance with (2.1) by setting 
^uu — Vfiv and the normalization of the kinetic term for Xs is standard. We shall compute the 
mass renormalization of Xs due to the one loop diagram of Fig. [5j More specifically we shall 
compute the contribution from the region of loop momentum integration where the momentum 
a, carried by the vector fields is small. In this limit we essentially compute the correction to 
the mass due to the Coulomb field of the x particle and the result is independent of the spin of 
X- Indeed it is easy to verify that our results remain unchanged if instead of x ^^ use another 
field with a different spin e.g. a Dirac field. This is important since the SO (32) spinor states 
for which we need to compute the mass renormalization in string theory are not scalars but 
transform in the 84+44 representation of the little group 5*0(9). 

In order to compute the contribution from the graph shown in Fig|5]we first need to choose 
a convention for the SO(32) generators T". In the vector representation we choose them to be 
32 X 32 matrices whose mn element is z, nm element is —i and whose other elements are zero. 
Collection of all such matrices for 1 < m < n < 32 form the basis of SO (32) generators. With 
this convention we have Try{T"'T^) = 26ab and hence the gauge kinetic term is normalized 
as —^Ff^j^F""^^ . As a result the gauge field propagator of momentum d in the Feynman gauge 
will be of the form —2i6abgfiu/i^^ ~ ^^)- The coupling of the gauge field to x in the Lgrangian 
density is of the form iA"-'^{T"-)ssi{x*s(^nXs' — d^x*sXs') from which we can read out the vertex. 
This gives the contribution to 5M'^ to be 

where the momentum integration is still written in the Lorentzian space, and {T°'T°') spinor 

21 



denotes the eigenvalue of T^-T"- on the spinor representation. With the normahzation of T" 
we have chosen, each T" has eigenvalue ±1/2 on the spinor representation and hence, after 
summing over the 496 generators of SO (32), T"'T"' will have eigenvalue 124 on the spinor 
representation. 



We shall analyze the contribution to (B.2) from the region of integration where i is small so 
that the low energy effective field theory description makes sense. In this case we can replace 
the vertex factor {2k^ — £^){2k^ — i'^) by Ak^k'^ = — 4(fco)^ and take it out of the integral. 
Making the Wick rotation £° — )• ii^ and denoting by is the Euclidean momentum we get 



In order to make connection with the string theory result (4.13) we now need to express the 



propagators in the Schwinger proper time formalism. We write 

'0 ^0 

Changing variable to T2 = 1/2 + -^2 and Z2 we get 

ilOn roo /•T2 



SM'„ = QQHk,rrr^j ^J dy2 j d., e-^^4-.((^.-^-m/^ (B.4) 



= 992{ko)'n^ I T-% /" ^^2 H dz2e~^^''^'^~^''/^'^"+-'"''"/^'-^''^''"+^'"^ 
J (271") Jo Jo 

^00 l'T2 

= 992(A;o)'7r2(27r)-i° / rfra / rf^2 (ts)"' e-^"^'/r2 ^ (35) 

Jo Jo 

where in the last step we have carried out the integration over Ie and have also used the 
on-shell condition k"^ = — M^ = —4 in the exponent. 



We shall now try to reproduce the same integral from the string theory result (4.13). We 



shall focus on the region of integration where both Z2 and T2 are large but Z2 remains small 



compared to T2. Since the integrand in (4.13) is invariant under z —^ r — z, there will be an 
identical contribution from the region where {t2 — Z2) is small compared to T2, and the effect of 
this will be to double the contribution from the ^2 << ^2 region^] Now since in computing the 
field theory contribution we examined only the graphs with bosonic intermediate states, we 
must do the same in string theory. This corresponds to restricting the sum over spin structure 

^Physically z — > r — z transformation exchanges the role of the two arms of the torus between the points 
and z. 

22 



V to 00 and 01 sectors only. We now use the following approximations to the various factors 
in (4.13) for large r2, ,22 y 



^^x{zf ~ 1-4 e^*^~2^*^ - 4 e^^^+2^^^ + ■ ■ ■ 



-16 



?9oo(^/2) -1 + 16 e~"*"+"^^ + 16 e"""""*^ + i20e^'"*"+'"*" + 



-16 



^01 (^/2) -1-16 e-""+"*^ - 16 e""'"-"*^ + i20e-'"*"+'"*" + 



16 



^10(^/2) ^ 



-16 



^n(z/2) ~ 

~ e"^"/^e-"" (1 + e"2"^" + e"^"^"" + ■ ■ ■ + e^'^^""^) + ■■■), 



{r]{T))-' ~ e"*"/^'(l + e-''^*" + ■■■). (B.6) 



Substituting these into (4.13), and taking into account the extra factor of 2 due to z ^- r — z 
symmetry, we get the string theory result for 6M'^ from the region where Z2 and T2 are large: 



5M'' c^ 8xAf{kof f Sri d^ze^^'^^-'^ 

Jt'2»1 Jl«Z2«T0, 



't2»1 Jl«Z2«T2 

1 + 4e~2"^" + ■ ■ ^(l + 4e-2"^"" + 4e-2-^-+2-i^ + . . .)((i + 1206^^"^"+^""" + ■ ■ ■) 
{1 + ■ ■ ■} {(1 + 2e~^^'~' + 2e-2'^^("-^) + ■ ■ ■) (1 + Me'^'^*" + ■ ■ .)e-^"^2/"2 (ra)"^ . 

(B.7) 

This expression has been organized as follows. After factoring out the leading terms inside each 



of the curly brackets in (4.13), we have written inside the three curly brackets the subleading 



terms from the terms inside the three curly brackets in (4.13). In particular since the leading 
holomorphic term has no factors of e"^'^*'^ or e"^'^*^ we have dropped all the subleading pieces 
containing factors of e^'^*'^ or e^'^*^. 



''Note that the leading term in the integrand has a stronger singularity in f compared to r and hence in the 
expansion in e^^'^*'^ and e~^'^'^ we must keep more terms than in their holomorphic counterparts. The general 
rule is that while expanding the holomorphic terms we only keep the leading terms except in those inside the 
sum over the spin structure v where we have to keep the first subleading terms. On the other hand for the 
anti- holomorphic terms we need to keep up to subleading terms of order e"^'^*'^, 
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Now for large T2 the integrals over Ti and Zi run from —1/2 to 1/2 without any restriction. 
Carrying out these integrals projects us into those terms which do not have any factors of 



^-2ttit qj, e"^'^*^. Thus we need to identify such terms in (B.7). Furthermore since we are only 
interested in picking up the contributions due to gauge boson exchange we need to consider 
only the subleading terms coming from the expansion of the terms inside the first curly bracket 
since this is what arose from the 32 left-moving fermions. The Q-'^'^^i'^-^) term in this expansion, 
which cancels the overall multiplicative factor of e^'^*^^"^), has coefficient (120 + 4) = 124. Thus 



the gauge boson exchange contribution to (B.7) is given by 



8 X 124Ar(A;o)^ f (fr /" t/^^e-^'^'Aa (^^)-5 ^ ^gg^ 



Comparing this with the field theory result (B.5) we now get 

992 



8 X 124 ^ ^ ^ ' 
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